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ABSTRACT 


It  is  shown  that  chemically  sensitive  field  effect  transistors  can  be 
used  to  study  stochastic  processes  at  the  solution/solid  interface  which  is 
not  tit  equilibrium.  Auto  spectral  densities  and  coherence  functions  for 
corrosion  of  aluminum>sillcon  alloy  at  xero-net  current  conditions  have  been 


obtained. 


INTRODUCTION 


Noise  in  chemically  sensitive  field  effect  transistors  (CHEMFET) 
originates  both  from  the  solid  state  part  of  the  device  and  from  the 
electrochemical  processes  at  each  interface.  In  our  previous  work  we  have 

so* 

shorn  that  the  analysis  of  noise  of  CHEMFET*  with  various  ion  selective 
ammbranes  at  chemical  equilibrium  yields  an  information  about  the  exchange 
currents  at  the  solution/membrane  interface  and  about  the  bulk  impedance  of 
the  membrane  itself  [l]. 

High  signal-to-noise  ratio,  small  input  capacitance,  small  surface  area 
of  the  electrode  and  low  parasitic  impedance  are  the  attractive  features  of 
CHEMFETs  for  the  study  of  stochastic  processes  in  electrochemical  systems. 
Because  each  chip  contains  at  least  two  transistors,  it  is  possible  to  use 
multichannel  analysis  methods.  These  methods  can  yield  an  unique  additional 
information  about  the  whole  system  and  they  reduce  contributions  from  the 
solid  state  part  of  the  system. 

In  the  first  part  of  this  series,  we  shall  outline  the  theory  of  noise 
analysis  in  CHEMFETs  and  describe  the  experimental  arrangement  which  was  used 
in  this  and  the  subsequent  studies. 


PRINCIPLE  OF  OPERATION 


The  basic  operation  of  CHEMFETs  has  been  described  elsewhere  [2].  The 
output,  drain-to-source  current  Ip  ,  for  a  n-channel  transistor  operated  in 
the  fion-saturated  region  [Vp  <  (VG  -  VT>]  is: 


h  =  K  f<VG  -  V 


(1) 


and  in  the  saturation  region  where  Vp  >  (V^  -  VT): 

JD  *  K'  (VG  “  VT)2 


(2) 


The  constants 
properties  of 
transistor  is 


K  and  K'  are  characteristic  of  the  geometry  and  material 
the  transistor  shown  in  Fig.  1.  The  transconductance  of  the 
defined  as 


g 


m 


(3) 


Because  the  magnitude  of  the  solid  state  part  of  the  noise  depends  on  the 
value  of  the  drain  current,  it  is  preferable  to  operate  the  transistor  in  so 
called  feedback  mode  in  which  the  average  drain  current  is  kept  constant.  In 
this  case,  the  gate  voltage  fluctuations  v  (t)  are  transformed  into  the  drain 


current  fluctuations  i(t)  through  the  relationship: 


i(t)  ■  g  v  (t) 
m 


(4> 


FLUCTUATION  ANALYSIS  THEORY 

In  this  section  a  brief  review  of  theory  relevant  to  the  present 
experiments  is  given.  For  a  fuller  account  of  the  theoretical  background  the 
reader  is  referred  to  textbooks  [3-5].  The  drain  current  in  the  channel  l^Ct) 
of  an  1SFET  can  be  expressed  as  [l]: 

ID  <t)  =  I  ♦  i(t)  (5) 

Where  I  is  the  mean  value  of  the  current  and  i(t)  is  the  fluctuating  part  of 
the  drain  current  (noise).  In  the  remaining  part  of  this  text  we  will 
consider  the  fluctuations  in  the  drain  current  only  and  for  the  sake  of 
notation  describe  drain  current  fluctuations  as  x(t)  or  y(t). 

When  one  is  interested  in  the  origin  of  the  fluctuations, 
characterization  of  the  signal  (by  means  of  parameters  for  instance)  is 
necessary.  Before  any  meaningful  signal  characterization  can  be  attempted, 
the  signal  has  to  be  qualified. 

This  qualification  can  be  carried  out  in  the  following  way:  Figure  2 
shows  graphically  the  concept  of  a  Probability  Density  Function  (PDF)  which  is 
defined  for  the  signal  x(t)  as: 


p(x) 


lim 

Ax 


Prob[x  <  x(t)  <  x  •»  Ax] 
Ax 
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In  general  p(x,t)  is  a  function  of  time.  If  thia  ia  not  the  cate,  that 
ia  when  p(x,tj)  *  pCx.t^)  for  all  tj,  tj,  it  ia  aaid  that  the  fluctuating 
signal  is  stationary.  By  its  definition  the  PDF  p(x,t)  is  normalised,  i.e., 
p(x,t)dx  ■  1,  where  the  integration  is  carried  out  over  all  values  of  the 
signal  x(t).  For  the  PDF  statistical  moments  can  be  calculated  in  the  usual 
way,  e.g.,  [6].  A  large  class  of  signals  is  formed  by  signals  with  a  Gaussian 
PDF  as  illustrated  in  Fig.  2.  When  a  fluctuating  signal  of  this  class  has 
'direct*  statistical  moments  defined  by: 

-mean  value: 


«  *<t>  >  -  |  J*  x<t>4t  ) 

-mean  value  of  the  square: 

<  x2(t)  >  =  i*™  (4  f  x2(t)dt  ) 


(6) 


(7) 


-variance : 

2  2  2 
o*  =  <  x  <t)  >  -  <  x(t)> 


(8) 


When  the  measuring  time  is  taken  long  enough  (usually  a  limit  T  — * 00  is 
Involved)  and  the  above  defined  moments  equal  the  corresponding  moments  of  the 
PDF  for  a  stationary  fluctuation,  the  signal  is  ergodic. 
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FREQUENCY  DOMAIN  ANALYSIS 


When  the  signal  is  stationary  and  ergodic  its  parameter*  are  independent 
of  the  observation  time.  This  property  allows  for  analysis  of  signal 
parameters  from  individual  time  records,  e.g.,  mean  and  variance  or  from  many 
different  records  as  e.g.  for  frequency  analysis  using  sampled  portions  of  the 
signal  of  finite  length  (to  be  referred  to  as  ensembles). 

Suppose  that  one  ensemble  is  represented  by  N  samples,  taken  at  time 
intervals  At;  N  is  usually  an  integer  power  of  2.  From  the  sequence  of 
samples  one  calculates  a  Discrete  Fourier  Transform  (DFT),  usually  using  a 
Fast  Fourier  Transform  algorithm  (FFT)  to  save  computation  time  [7],  From  the 
Fourier  components  estimates  of  spectral  densities  Gx(k)  are  calculated  as 
discrete  frequencies  kfQ,  unique  up  to  k  -  N/2  -  1.  The  sequence  G^Ck)  is  a 
discrete  representation  of  the  auto  spectral  density  function  for  the  signal 
X(t).  The  true,  continuous,  single-continuous,  single-sided  auto  spectral 
density  G^Cf)  is  defined  as: 

T 

Gx(f)  *  2  —  f  Xk  (f,  T)  2dT  (9) 

T  o' 

The  random  part  of  the  error  in  the  spectral  estimates  is  given  by 
Bendat  and  Piersol  [3]  for  the  case  of  a  signal  with  a  Gaussian  PDF  as  a 
normalized  standard  error: 
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(10) 


where  H  equals  the  number  of  degrees  of  freedom  in  statistical  sense  (each 
time  record  has  two  degrees  of  freedom  due  to  the  fact  that  the  spectral 
estimate  has  a  real  and  an  imaginary  component).  This  implies  that  the 
standard  deviation  in  a  spectral  estimate  calculated  from  one  time  record  is 
100%,  independent  of  the  length  of  the  record. 

The  number  of  degrees  of  freedom  can  be  increased  in  two  ways:  by 
smoothing  of  the  estimates  at  the  expense  of  spectral  resolution  or  by  segment 
averaging.  The  latter  technique  requires  a  longer  observation  time,  but 
preserves  the  spectral  resolution,  and  is,  therefore,  adopted  in  this  paper. 

The  bias  part  of  the  error  in  the  spectral  estimates  is  a  function  of 
the  ’shape*  of  the  spectral  density  function.  For  flat  spectra,  1/f  noise  and 
Lorentzian  spectra  this  error  is  neglibly  small  at  the  frequencies  of  interest 

[3]. 

From  the  definition  of  the  spectral  density  it  follows  immediately  that 
if  a  signal  can  be  considered  to  be  the  sum  of  some  stationary,  ergodic 
signals  the  spectral  density  of  the  sum  equals  the  sum  of  the  spectral 
densities  of  the  individual  signals: 


x(t ) 

G  (f) 
x 


i  =  1,  2,  3  .  .  . 
i  =  1,  2,  3  .  .  . 


(31) 


In  a  previous  paper  [l]  this  property  was  used  by  Haemmerli  et  al.  Co 
discriminate  between  noise  originating  in  different  parts  of  the 
ISFET/electrochemical  cell  system.  It  is  clear  from  eq.  (11)  that  as  soon  as 
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one  of  the  different  nolle  sources  is  stronger  than  the  others  this 
discrimination  becomes  impossible.  The  same  conclusion  was  reached  by 
Haenmerli. 

TWO  CHANNEL  APPROACH 

In  order  to  reduce  the  device  noise  in  a  ISFET/electrochemical  system,  a 
two  channel  approach  was  developed  based  on  the  following  reasoning:  Suppose 
that  two  ISFETs  have  an  electrochemical  noise  source  in  common  but  that  their 
'solid  state*  noise  is  not  correlated.  For  the  fluctuations  in  the  two  drain 
currents  we  can  write: 


x(t)  =  dj(t)  +  e(t) 

y(t)  =  d2(t)  -  e(t)  (12) 

dj(t)  is  the  device  noise  for  FET1,  d2(t)  is  the  device  noise  for  FET2  and 
e(t)  is  the  (electrochemical )  noise  common  to  both  channels.  The  true 
continuous  single-side  spectral  density  functions  are: 


Gx(f)  =  tDl(f)  *  >3[D1(f )  ♦  E(f )]* 

Gy(f)  =  [D2(f )  ♦  E(f)][D2(f)  +  E( f  )]* 

Gxy(f)  =  [D2( f )  4  E(f JJtDjCf )  4  E( f ) ]* 


(13) 
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fi¬ 

ll 


£ 


r 


where  (*)  denotes  a  complex  conjugate.  D^(f)  and  l>2(f)  are  Fourier  transforms 
of  signals  dj(t)  and  d2(t),  E(f)  is  the  Fourier  transform  of  the  noise  e(t) 
common  to  both  channels.  If  we  drop  the  frequency  dependence  for  convenience 
of  notation,  expanded  eq.  (13)  becomes: 


G 

x 

G 

y 


*  *  *  * 

D}D1  4  EDj  4  D^E  4  EE 

■*  *  *  * 

D2D2  4  ED2  4  D2E  4  EE 

*  ★  it  ic 

D2D1  4  D2E  4  EDj  4  EE 


(14) 


If 


d  ,  d2  and  e  are  uncorrelated  noise  sources  all  cross-terms  in  eq.  '14): 


EDj*  =  DjE*  =  ED2*  =  D2E*  =  DjD2*  = 


W 


are  rero  and  hence: 


*  ★ 

DjD  4  EE 

*  * 

D2D2  4  EE 


EE 


(15) 


The  cross  spectral  density  function  is  thus  a  function  of  electrochemical 
noise  common  to  both  channels  only. 

A  derived  quantity  is  the  coherence  function  defined  as: 


yd) 


2 


=  G  (f)  C  (f) 


x  y 


(16) 


I 


Theoretically,  this  function  has  a  range  from  zero  for  no  coherence  to  one  for 
complete  coherence.  It  can  be  viewed  as  a  measure  for  the  amount  of  signal 
common  to  both  channels. 

The  same  considerations  as  above  apply  to  the  statistical  errors  in 
estimates  of  the  auto  spectral  density  functions  obtained  by  Fourier  analysis 
of  time  records.  Error  estimates  in  the  cross  power  spectral  density  and  in 
the  coherence  function  are  dependent  on  the  amount  of  coherence  between  the 
two  signals  and  are  relatively  complicated  [3,8].  The  number  of  degrees  of 
freedom  is  important  especially  at  low  levels  of  coherence  as  Table  1  shows: 


TABLE  1 


Error  estimate  in  coherence  function  at  99 X  confidence 
level  and  for  coherence  »  0.01 


Degrees  of  Freedom 

Lower  Limit 

Upper  Limit 

100 

0 

.06 

1000 

.002 

.022 

10000 

.007 

.013 
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EXPERIMENTAL 


The  transistor  chips  used  in  this  study  are  identical  to  those  used 

previously  [l].  They  were  mounted  on  1.5  mm  diameter  PVC  tubing  and 

encapsulated  with  a  high-grade  epoxy  resin  (Epon  826  with  Jeffamine  D230 

crosslinker).  The  two  gates  used  for  the  cross-correlation  measurements  are 

connected  as  shown  in  Fig.  1.  In  this  study  we  have  used  25  pm  diameter  Al-Si 

(1%)  wire,  the  geometrical  area  of  the  metal  exposed  to  the  solution  was 

—6  2 

approximately  7  x  10  cm  .  During  etching  experiments  this  area  was  renewed 
by  cutting  the  sample  perpendicularly  to  the  wire  axis  after  each  experiment. 
Unless  stated  otherwise,  the  measurements  were  done  at  room  temperature.  The 
actual  temperature  at  the  sensor  site  was  determined  by  pre-calibration  of  one 
of  the  MOSFETS  available  on  the  transistor  chip. 

The  noise  measurements  were  done  under  two  sets  of  experimental 
conditions:  in  pH  9.0,  0.1  M  borate  buffer  (no  corrosion)  and  in  etching 
solution  consisting  of  1600  ml  H^PO^  (35%  Fisher)  100  ml  HNO^  (697.  Fisher), 

100  ml  acetic  acid  (glacial,  Fisher)  and  200  ml  of  water.  This  solution  was 
diluted  with  distilled  deionized  water  for  concentration  dependence  studies  as 
necessary.  The  solutions  were  degassed  with  nitrogen  and  stored  in  glass 
bottles  under  nitrogen.  They  were  delivered  to  the  transistor  under  constant 
pressure  using  the  manifold  shown  in  Fig.  3.  The  volumetric  flow  rate  was 
controlled  by  adjusting  the  nitrogen  pressure.  Because  of  the  irregular 
geometry  of  the  encapsulated  sample,  the  absolute  face  velocity  of  the  flowing 
solution  at  the  sample  surface  could  not  be  determined.  The  electrical 
circuits  used  for  monitoring  the  drain-to-source  current  are  shown  for  a  two- 
channel  operation  in  Fig.  4.b.  The  output  signals  were  amplified  with  low 
noise,  variable  gain  (1-10  K)  pre-amplifiers  (Princeton  Applied  Research, 

Model  PARC113A)  and  analyzed  using  Hewlett-Packard  Signal  Analyzer  (HP  3582A). 
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The  accessible  frequencies  of  chis  measuring  setup  were  between  0.8  Hz  and  25 
kHz. 

RF.SULTS 

In  the  preliminary  experiments  we  have  used  diluted  hydrochloric  acid 
(0.1  mM  to  0.1M)  as  the  reaction  medium.  For  concentrations  above  5mM  HC1  the 
reaction  was  too  violent  yielding  bubbles  of  hydrogen  which  created  problems 
in  the  flow  cell.  For  concentrations  below  1  mM  the  change  of  surface  pH  due 
to  low  buffering  capacity  led  to  irreproducible  results.  Nevertheless,  the 
auto  spectral  densities  in  these  experiments  exhibited  a  pronounced  maximum  in 
the  region  between  1-10  Hz  which  increased  with  increasing  concentrat ion 

—1 8  1  —1  i 

of  the  acid  from  5  x  10  A  Hz  for  0.1  mM  HC1  to  5  x  10~  A  Hz~  for 
0.1  It  HC1 

In  order  to  control  the  surface  concentrations  better  a  buffered 
solution  normally  used  for  etching  aluminum  in  the  integrated  circuit 
fabrication  process  [9]  has  been  used  in  various  dilutions.  For  undiluted 
solutions  at  room  temperature  the  typical  etching  rate  is  less  than  1000  K 
min  .  Therefore,  for  the  duration  of  our  experiment  (~10  min.)  the  removal 
of  the  material  has  not  significantly  altered  the  flow  pattern  around  the 
sample . 

The  auto  spectral  densities  for  etching  solutions  of  various  dilutions 
are  shown  in  Fig.  5.  These  are  the  results  of  123  averages.  There  is  a  clear 
evidence  of  a  Lorentzian  character  of  the  spectrum  (slope  -2)  for  the  two 
highest  concentrations.  The  intensity  of  the  spectrum  increases  with  the 
increasing  concentration  of  the  etching  solution  in  general  agreement  with  the 
results  obtained  with  hydrochloric  acid. 

The  coherence  function  (eq.  16)  for  the  data  presented  in  Fig.  5  is 
shown  in  Fig.  6.  It  again  shows  the  expected  concentration  dependence.  In 
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buffer  solution  where  no  significant  corrosion  takes  place,  its  value  is  0.025 
and  Independent  of  frequency  up  to  1  kHz.  The  fact  that  the  coherence  has  a 
non-zero  value  indicates  the  presence  of  common  electrochemical  noise  which 
probably  originates  from  the  pseudo-exchange  current  (non-equilibrium)  which 
takes  place  at  the  metal/solution  Interface.  This  statement  has  to  be  taken 
within  the  context  of  the  statistical  accurracy  of  the  coherence  function 
(Table  1). 
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CONCLUSION 


In  this  paper  we  have  outlined  the  framework  for  study  of  non- 
equilibrium  fluctuation  phenomena  in  chemically  sensitive  field  effect 
transistors  both  from  the  theoretical  and  from  the  experimental  point  of  view. 
The  corrosion  of  aluminum/silicon  alloy  is  used  only  as  an  illustration  that 
this  approach  is  feasible  and  no  interpretation  of  the  experiments  is 
attempted.  The  corrosion  of  metals  and  semiconductors  is  commonly  studied  by 
conventional  electrochemical  means  in  which  the  important  corrosion  parameters 
are  derived  from  current-voltage  characteristics  which  are  extrapolated  to 
aero  net  current  conditions.  Bertocci  [10-12]  has  introduced  the  concept  of 
coherence  to  corrosion  measurements.  It  has  allowed  him  to  determine  the 
causal  relationship  between  current  and  voltage.  The  fundamental  difference 
in  our  measurements  is  that  they  are  performed  at  zero-net  current  condition. 
Thus,  the  only  independent  variables  are  surface  concentrations,  temperature 
and  pressure.  This  approach  is  possible  because  the  coupling  between  the 
electrochemical  process  and  the  Interface  and  the  observable  (fluctuation  of 
the  drain-to-source  current)  is  through  the  electric  field.  In  the  following 
paper  we  will  show  that  this  mode  of  operation  yields  an  unique  approach  to 
the  study  of  electrochemical  reactions  of  various  materials,  including 
insulators . 
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Figure  1.  Schematic  diagram  of  the  transistor  chip  (drawing  not  to  scale). 

D  -  drain,  S  -  source,  I  -  insulator,  M  -  metal,  E  -  encapsulate, 

R  -  reference  electrode,  and  SB  -  substrate.  is  drain-to-source 
current  and  is  the  applied  gate  voltage 

Figure  2.  Illustration  of  Probability  Density  Function  (PDF)  for  an  arbitrary 
signal,  (a)  Time  record,  (b)  PDF. 

Figure  3.  Experimental  arrangement.  A,  the  amplifier,  P,  battery  pack,  T, 
the  transistor  under  study,  REF,  the  reference  electrode,  F,  the 
flow-cell,  H,  the  heater,  S,  the  solution  inlet,  W,  the  outlet  to 
waste,  0,  the  amplifier  output,  V's  are  control  valves,  L,  the 
hydrodynamic  load,  E,  the  etching  solution,  and  B,  the  buffer 
( - ),  Faraday  shielding. 

Figure  Circuit  used  in  two-channel  measurements;  REF  -  reference 

electrode;  E  -output  voltage;  D  -  drain;  S  -  source. 

UU  T 

Figure  5.  Autospectral  density  for  corrosion  of  aluminum  at  three  different 

concentrations  of  etching  solution.  Vertical  axis  current  density 

2 

units  are  in  A  /Hz. 

Figure  6.  Coherence  function  obtained  in  two  channel  mode  for  three  different 
concentrations  of  etching  solution  and  for  buffer  solution. 
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SOLUTION 


LOG  AUTO  SPECTRAL  DENSITY 


LIST  OF  SYMBOLS 


o(f) 

E(f) 

d(t) 

«U) 

v_ 


D 

i(t) 

K,K' 

v(t) 

x(t) 

y(t) 

p(x,t) 

<x> 

<x2> 


Gx(f) 

Gy(f) 

g  a> 
*y 

N 

M 

Y2(f ) 


Fourier  transform  of  device  noise 
Fourier  transform  of  electrochemical  noise 
Device  noise  in  time  domain 
Electrochemical  noise  in  time  domain 
Drain-source  voltage 
Gate-source  voltage 
Threshold  voltage 
Drain  current 

Average  value  of  drain  current 
Fluctuation  in  drain  current  in  time  domain 
Transconductance  of  FET 
Constants  of  FET 
Fluctuations  in  gate  voltage 
Fluctuating  signal  in  time  domain 
Fluctuating  signal  in  time  domain 
Probability  density  function  for  x(t) 

Mean  value  of  signal  x(t) 

Mean  value  of  the  square  for  x(t) 

Variance  of  x(t) 

Auto  spectral  density  for  signal  x(t) 

Auto  spectral  density  for  signal  y(t) 

Cross  spectral  density  for  signals  x(t)  and  y(t) 

Number  of  samples  in  finite  length  discrete  time  record 
Number  of  averages  in  spectral  density  estimate 
Coherence  function 
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